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Комплексний контроль техначного стану 
суднового двигуна з турбонаддувом 


В статье анализируются структура и характеристики разрабатываемого модуля удаленного контроля 
частотных характеристик судового двигателя с газотурбинным наддувом (СДВС). Модуль предназначен 
для непрерывного дистанционного контроля частоты вращения коленчатого вала (КВ) и частоты 
вращения ротора газотурбонагнетателя (ГТН), а также комплексных характеристик технического состояния 
СДВС: разности цилиндровых мощностей и величины остаточного небаланса вибрации ротора ГТН. 
Принцип работы модуля основан на анализе спектральных характеристик виброакустических сигналов 
комбинированного двигателя. В модуле реализуется процедура БПФ и численно решается задача 
устранения эффекта «утечки», возникающего в результате дискретного преобразования Фурье. 
Ключевые слова: судовой дизель, диагностика, спектральный анализ. 
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Введение 


С внедрением на флоте технологий интернет-мониторинга технического состоя- 
ния энергетических установок стал актуальным вопрос непрерывного постоянного 
контроля их основных параметров независимо от штатных приборов. Большинство 
главных судовых дизелей (СДВС) оборудованы тахометрами для определения частот 
вращения коленчатого вала и ротора ГТН. Датчики тахометров непосредственно свя- 
заны с коленчатым или распределительным валом двигателя и валом турбины наддува. 
Принцип действия датчиков различен — оптические, индуктивные, датчики Холла и др. 
Они осуществляют непосредственное измерение частот вращения указанных объектов. 
Сигналы от тахометров выведены на панель приборов центрального поста управления и 
в систему аварийно-предупредительной сигнализации. Поскольку эти частотные парамет- 
ры являются основными характеристиками эксплуатируемого двигателя, всю схему из- 
мерения заверяет выбранное судовладельцем сертификационное общество (Регистр) [1]. 
Использовать гальваническое присоединение в регистровую схему измерения для пере- 
дачи информации во внешние контролирующие устройства в большинстве случаев 
запрещено. Поэтому в качестве альтернативы непосредственным измерениям пред- 
лагается использовать виброакустический способ определения частот вращения КВ 
и ГТН с помощью методов спектрального анализа. Кроме частотных характеристик 
в спектрах виброакустических сигналов содержится диагностическая информация, 
позволяющая производить комплексный анализ технического состояния СДВС и ГТН. 


Виброакустический метод определения частот вращения 
двигателя и газотурбонагнетателя СДВС 


На судовых дизельгенераторах в отличие от главных дизелей в большинстве слу- 
чаев не устанавливаются тахометры турбин наддува и виброакустический контроль 
частоты их вращения является удобным для практического применения методом. 
С помощью методов спектрального анализа на базе быстрого преобразования Фурье 
(БПФ) и последующего численного решения задачи устранения «утечки» в дискретном 
спектре можно снизить погрешность измерения частотных характеристик комбиниро- 
ванного судового дизеля. Главной целью применения такого метода является непре- 
рывный дистанционный контроль частот КВ и ГТН в задаче опПпе-мониторинга ре- 
жимов работы судовых дизелей и возможность проверки показаний штатных приборов. 
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Рисунок 1 — Спектр звукового сигнала компресора ГТН главного дизеля 6180МСЕ 
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В [2] было показано, как с помощью амплитудного спектра виброакустического 
сигнала ГТН определить частоту вращения ротора турбины. Анализ литературы [1-5] 
показывает, что в спектре виброакустического сигнала ГТН максимальная гармоника 
находится на лопаточной частоте воздушного компрессора (частота вращения ротора 
умножить на количество лопаток компрессора, рис. 1). 

Анализ рис. 1 показывает, что необходимо вводить ограничения на поиск макси- 
мальной гармоники, соответствующей лопаточной частоте ГТН. Для этого необходимо 
знать частоту вращения ротора ГТН на номинальном режиме и количество лопаток 
компрессора. Так, для малооборотного главного судового дизеля (МОД) 6180МСЕ 
частота вращения ротора турбины на номинальном режиме равна 9000 грт. Тогда 
ожидаемая частота лопаточной гармоники на номинальном режиме: 


Еа(к) = 9000 грит / 60 * 20 = 3000 Ги. 


Минимально и максимально возможные частоты вращения ГТН и, соответст- 
венно, Частоты генерируемых сигналов можно принять следующими: 


ти -> 6000 грт / 60 * 20 = 2000 Гц; тах -> 15000 гри / 60 * 20 = 5000 Гц; 


На рис. 1 видно, что в указанном диапазоне существует только один четко выра- 
женный максимум. Зафиксированная по спектру сигнала частота вращения ротора 
ГТН (в рамках принятых ограничений) составила: 


ТОК тра = 2948 Гц / 20 * 60 = 8844 грт. 


Анализ литературы [1-5] и проведенные натурные испытания показывают, что 
в спектре виброакустического сигнала двигателя, в области низких частот (до 100 Гц) 
выделяется гармоника, соответствующая цилиндровой частоте вибрации корпуса. 
Для двухтактных дизелей МОЛ: 


Ед су1(2) = грт / 60 * 1су! 


Для четырехтактных дизелей СОД и ВОД, с учетом одного оборота КВ на рабо- 
чий цикл и одного оборота на продувку цилиндра: 


Ра су1(4) = (грт / 60 *1су0 /2 


В связи с дискретным представлением сигнала в спектре появляется несколько 
кратных гармоник и в этом случае обязательно необходимо вводить ограничения 
«сверху» и «снизу» на поиск цилиндровой гармоники. Алгоритм формирования огра- 
ничений очевиден: вычисляется в герцах минимальная и максимальная возможные 
цилиндровые частоты, исходя из частоты вращения КВ на номинальном режиме. 
Так, для дизельгенератора 5и7ег 8АГ.25/30, работающего по нагрузочной характе- 
ристике, частота вращения поддерживается равной 750 грт. 


Еа су = (750 / 60 * 8) /2 = 50 Гц (рис. 2) 


С учетом максимально возможного кратковременного отклонения частоты не 
более +5% на рис. 2 показаны ограничения для поиска цилиндровой гармоники. 

Испытания дизельгенератора были проведены на 4 нагрузочных режимах: 240, 
280, 320 и 350 кВт. В спектре сигнала были определены цилиндровые и лопаточные 
гармоники, с последующим уточнением их амплитуды и частоты путем устранения 
эффекта «утечки» в дискретном преобразовании Фурье (ДПФ). В результате были 
построены графики зависимости частот вращения КВ и ГТН от нагрузки в исследо- 
ванном диапазоне (табл. 1, рис. 3, 1 су = 8, 1 лопаток ГТН = 20). 
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ВРМ(кв) = 60 * Еа(ц)*2 /тсу| (4-тактн.) 
ВРМ(кв) = 60 * Ра(ц) /1су1 (2-тактн.) 
ВКРМ(кв) = 60 * Ра(л) / 1 лопаток 


Таблица 1 — Испытания дизельгенератора 


№ кВт | Еа(@® КРМ Еа(л) КРМ 
(КВ) (ГТН) 

240 ВИ 175.5 3508 10524 
280 51,6 774 7 51 
320 51,6 774 3950 11850 
350 51,5 12 4045 12135 
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Рисунок 2 — Спектр звукового сигнала дизельгенератора СОД 8АТ25/30 
на нагрузочных режимах 240, 280, 320, 350 кВт 
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Рисунок 3 — Мощность 8АТ.25/30 и частота вращения КВ 
в зависимости от частоты вращения ГТН 


«Штучний 1нтелект» 32012 369 


6в Варбанец Р.А., Ивановский В.Г., Кучеренко Ю.Н., Головань А.И. 


Необходимо отметить, что статическая характеристика регулятора частоты вра- 
щения дизеля 8АТ.25/30 имеет наклон, который можно численно оценить. Это видно 
по уменьшающейся частоте вращения при увеличении нагрузки генератора. Зависи- 
мость мощности от частоты вращения ГТН, в выбранном диапазоне нагрузок, слабо 
квадратичная и, можно сказать, почти линейная. 


Устранение эффекта «утечки» 


В процессе анализа дискретного спектра виброакустических сигналов для оценки 
их частотных и амплитудных характеристик приходится решать задачу устранения 
эффекта «утечки». Этот эффект является следствием конечности анализируемой вре- 
менной реализации и ее дискретного представления через АЦП. В качестве примера на 
рис. 4 представлены амплитудные спектры одного и того же синусоидального сигнала 
(А = 10000 ед.) с целым (а) и нецелым (6) числом отсчетов на один период сигнала. 

Эффект «утечки» или вытекания мощности из спектральных пиков на соседние 
спектральные линии (рис. 4) считается одной из главных погрешностей дискретного 
преобразования Фурье. Если частоту сигнала представить в виде 

у = М/Т, 

где Т — период сигнала, а М =п + о, где п — целое, а 0 < ос < 1, то максимальные 
искажения амплитуды, частоты и фазы центральной гармоники и утечка мощности в 
соседние будет наблюдаться при о = 0,5 [6]. 

Таким образом, если анализировать параметры исходного сигнала по его спектру, 
т.е. по центральной гармонике, то полученные амплитуда, частота и фаза будут иска- 
жены в случае нецелого количества отсчетов сигнала на его период. На практике для 
дискретной записи сигналов используют АЦП с выбранной и зафиксированной час- 
тотой дискретизации. Понятно, что число отсчетов на период никогда не будет целым, 
а значение с будет меняться от 0 до 1 в зависимости от собственной частоты измеряе- 
мого сигнала, при этом будет меняться точность оценки параметров сигнала по цент- 
ральной гармонике. 
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Рисунок 4 — Эффект утечки ДПФ 
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Для уменьшения эффекта «утечки» наиболее широкое распространение получил 
метод оконных преобразований. Суть метода проста: для уменьшения разрывов на 
краях ряда с целью ослабления утечки нужно уменьшить амплитуду сигнала возле 
краев. Такое масштабирование осуществляется в ходе умножения реализации на окно 


специальной формы х,„’ =х, М (Г, например: 
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Таблица 2 — Оконные функции, применяемые для уменьшения эффекта «утечки» 


и 
Окно Хемминга (п) = 0.53836 — 0.46164 сс ( м .) 


2 
Окно Ханнинга (п) = 0.5 (1 соя ( ыы )) 


М-1 


Окно Кайзера 


В результате применения оконных функций спектр исходного сигнала может 
сильно исказиться, но при этом уменьшится зависимость амплитуды, частоты и фазы 


от значения о. Это значит, что центральную гармонику х. можно использовать для 


приблизительной оценки параметров сигнала с некоторой постоянной погрешностью, 
которую можно учесть. Точнее можно устранить эффект «утечки» численным методом, 
основанным на обработке комплексных результатов ДПФ. В [6] предложено уточнить 
частоту т, фазу ф и амплитуду А исходного сигнала по значениям двух максималь- 


ных гармоник в спектре Х 4 к: Для этого предлагается численно решить систему 
комплексных уравнений: 


|2 (т‚ф), 12 (т,Ф) и — |х, Ха 
Ат8(Ё(т,ф),) = Атз(Х,) 


где параметры (-й гармоники: 
Х, = Ке, + 7 Пт, ; 


Х, = МА, е**, 
А, = су Вел Пт?» , 


[п 
Фа =атав (> 


(0 


=) = А’в(Х,). 
о, 

Коэффициенты гармоник можно представить в виде: Х, = (А, /2)7(т,ф), ‚ где 
7(т,Ф), — комплексная функция, не зависящая от амплитуды, но зависящая от час- 
тоты и фазы: 

В ве о-в —1 й вчтея —1 
(т), 262’ едю 1 ыбею — 


е м -1 е м“ -1 
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Систему уравнений (1) нужно решать в том случае, если гармоники слева и 
справа от центральной не равны нулю (больше заданной малой величины д, рис. 46): 


ОАО 


Если Х,,=0,Х,., =0, то эффект утечки отсутствует, и частота, амплитуда и фаза 
центральной гармоники соответствуют параметром измеряемого исходного сигнала. 

Авторами была решена представленная выше система двух комплексных урав- 
нений (1) и написана программа, демонстрирующая решение проблемы. В первую 
очередь, необходимо было записать выражения для действительной и мнимой частей 
комплексной функции 7(т,Ф), : 


уу Уи 
Вели) = г + В. 
1 2 


Уи И+Е 


п — 
2(т,ф) р) р) 
1 2 


) 


где 


о пт 2^\ \ элт 2 
р, =| со чт += |, ВБ, =| с09 += |-1| + ——+—— |, 
мм мм мм мм 


скофорко аки) эти} [оо и | 


= (- зи (со{—=2лт-+=2^)-- Оенсофо или ВР И 
= (со {со$2лт- 2л^)-—1)— и (ф)за( (2пт-+ 2лк)) › ыы и 
) 


= (и (со{2лт-- 2) 1) +со{ф)зи(2лт-+ 2лк)) хз 


= (со$фсо<— 2т-+ 21) 1)+ зтФ )51(- 2жт-+ 2л^)) х ние к 
У = (-1®Ф\со{2лт+ 2) -1)- со{ф)зщ2лт-+2лк)) < ик й 


= (со$ф(со{2лт+ 2лк) в |- зш(ф)зт(лт-+ 2лк)) х51 


Аргумент комплексной функции И(т,ф), : 


[п 
Ате(7(т,ф)) = поить 


(т.о) 


Модуль комплексной функции 7(т,Ф), : 


|2(т.Ф))| = био)? + Веди 
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Система (1) решается методом итераций. Вначале задается значение ф = Аге(Х,) 


и численно решается первое уравнение системы при одном неизвестном то. Получен- 
ное значение т задается во второе уравнение и оно решается относительно ф. Теперь 


полученное значение ф снова задается в первое уравнение, для получения уточненного 


значения т, ит.д. Признаком окончания решения будет незначительное на очередной 
итерации изменение тиф: 


т И (= 
у -Ф/(& | 
После завершения итераций амплитуда сигнала определяется из соотношения: 
А, = АБ5(2Х 12 (т,Ф),). 


Несмотря на кажущуюся громоздкость системы (1) и необходимость численного 
итерационного решения составляющих ее уравнений, необходимо отметить, что реша- 
ется она быстро. Для всех исследованных случаев требовалось не более 5 полных ите- 
раций для обеспечения заданной погрешности менее 0,5% по частоте и фазе. На рис. 5 
показано решение системы для ситуации сильно выраженного эффекта утечки о - 0,5. 
В данном случае потребовалось четыре полных итерации системы. Амплитуда и час- 
тота восстановлены до заданного в исходном сигнале значения (А1 = 0,8, \1 = 5,0) с 
точностью до 5-го знака после запятой. При этом амплитуда центральной (макси- 
мальной) гармоники в спектре после ДПФ до процедуры восстановления была равна 
А1’ = 0,521 (ошибка 35%!). Ошибка в оценке частоты исходного сигнала по частоте 
центральной гармоники тоже может быть значительной. Она зависит от частоты АЦП 
и частоты исходного сигнала. С увеличением частоты АЦП она будет уменьшаться. 


. тапоти 0000 >| А . ыИ 0,000 $ 10 |271 0,00 +| | пос МАХ: 0,521 
геад Не ро В 
А1 д2 .000 > 2 Е Е 0,000 > Г 22 0,00 < 384 $ 
А2 * зи2пиТ+Е] +2 | | = ги = ЕРО: 104,000 
АЗ '3:_ < 0400 000 м >< | ГЕО: 104438 


. . | Г зес/Н2 м 
00134А20040558413 ы © =1Е- 


К, К+1 = 104, 105 


^ 
3000 + ВеХК, К+1 = 0,521, -0,497 
|иХК, К+1 = 0,017, -0,016 


А = 0,032 

11=20,еди=-2.62Е-6, т =104,4875193,е 
ВЕЗЕТ =22.еци=-7.04Е-7, |=-85.917.е10 
11=20,еди=2.62Е-6, п =104,4898501 „ег 
1и=22 еди=-6,75Е-7, Р=-86,336,еи0 


1и=20,едч=-2,9Е-6, п =104,4898443.ет 
1и=22еди=-6,75Е-7, Р=-86,335,еп0 


и=20.еди=-2.35Е -6, т =104,4898443.е 
1=22.еди=-6.75Е-7, Р=-86.335,еи0 


'АМРЕ = 0.80000 

= = —= ‘1 = 5.000 [104,490] 
Ве0172 е0.52 2-04 РИ = 86.335 

_— 0172 т0.005 032 0521 т 0.017 0032 


Везий РЕТ 


3 *| 0144 Нг 
1023 >| 48,952 На 
#СН0) Г. 


А>0 м 
зака) м 
1" М 
2= 25.638 


Рисунок 5 — Окно программы демонстрации ДПФ и решения проблемы «утечки» 
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Решение системы (1) не связано с выделением дополнительной памяти для хра- 
нения объемных массивов данных и расчетных коэффициентов, как в случае с БПФ. 
В связи с этим алгоритм может быть запрограммирован на современном ОЗР конт- 
роллере, реализующем БИФ. Несмотря на итерационное численное решение (1), такая 
процедура восстановления незначительно увеличит общее время расчета и позволит 
получить не только спектр сигнала, но и восстановленное значение основной частоты 
и амплитуды исследуемого сигнала. 


Общие характеристики модуля 
виброакустического контроля параметров 
судового двигателя с турбонаддувом 


Модуль решает задачи определения частоты вращения КВ и ГТН двигателя по 
восстановленным параметрам дискретного спектра виброакустическоих сигналов дви- 
гателя и компрессора ГТН. Модуль работает на базе современного ОЗР контроллера, 
который реализует БПФ на базе 2048 точек и устраняет эффект «утечки» ДПФ. В ре- 
зультате кроме частот КВ и ГТН на выходе модуля производится непрерывный контроль 
амплитуд виброускорений двигателя на цилиндровой частоте и амплитуд виброуско- 
рений ГТН на основной частоте вращения. Первые амплитуды характеризуют разбаланс 
цилиндровых мощностей, а вторые отвечают за общую вибрацию ГТН на основной 
частоте вращения. Контроль уровня вибрации ГТН на частоте вращения — это комп- 
лексный контроль технического состояния ГТН, его подшипников и проточной части 
по газовой и воздушной стороне [7]. 

Модуль устанавливается на двигателе в районе ГТН. Настроечные данные и рас- 
считанные характеристики передаются по последовательному интерфейсу. 

Входные сигналы: 

1) виброскорость/ускорение блока цилиндров (до 100 Гц, датчик № 1); 

2) виброакустический сигнал компрессора ГТН (2 — 20 КГц, датчик № 2). 

Выходные параметры: 

1) частота вращения КВ двигателя (гртп); 

2) амплитуда цилиндровой гармоники (Асу|); 

3) частота вращения ГТН (ТОЕгрт); 

4) амплитуда гармоники на основной частоте вращения ГТН (Атов); 

Еа(ГТН) = Еа(лопаточная) / 1 лопаток 

Настроечные параметры: 

[) частота опроса датчика № 1 (1); 

2) минимальная КРМ двигателя (грип); 

3) максимальная ВРМ двигателя (грий»»х); 

4) тактность («0» -— 2-тактный, «1» — 4-тактный); 

5) количество цилиндров (1су)); 

6) частота опроса датчика № 2 (у2а‹); 

7) минимальная ВРМ ГТН (ТОВ трте»); 

8) максимальная ВРМ ГТН (ТУВгрти»»); 

9) количество лопаток компрессор ГТН (1); 

10) тип применяемых в ДПФ окон («0» — без окон, «1» — Хемминг, «2» — Блэкмен, 
«3» — Кайзер); 
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11) метод устранения эффекта утечки («0» — трехточечная аппроксимация [5], 
«1» — итерационный алгоритм); 

12) адреса каналов двигателя и ГТН (сБе», — первый полубайт, сти — второй 
полубайт). 


Заключение 


Виброакустический метод определения частотных характеристик судового дви- 
гателя с турбонаддувом имеет ряд преимуществ перед традиционными, предполага- 
ющими применение тахометров. В первую очередь это стоимость и простой способ 
инсталляции. Кроме того, постоянный анализ спектра виброскорости/ускорения блока 
цилиндров двигателя и ГТН дает возможность анализировать техническое состояние 
двигателя. Точность определения частот КВ и ГТН, по предварительной оценке, выше 
точности применяемых на практике тахометров, особенно при определении частоты 
вращения ГТН. Это объясняется тем, что происходит измерение не самих частот 
вращения, а частот, кратных количеству цилиндров или лопаток компрессора ГТН. 
Однако вопрос численной оценки погрешностей предлагаемого метода требует от- 
дельного решения. 

Виброакустический модуль определения частотных характеристик не является 
прямой альтернативой применяемых на практике тахометров, но может быть исполь- 
зован в составе системы непрерывного мониторинга и диагностики судового ДВС. 
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ВЕЗОМЕ 

В.А. УагБапеб, У.СЛуапоузКу, У.М. КисйегепКо, А.[. Со]оуап 
Моппоппе о{Магте ТитБо Г\езе! Епоте Тесйтса[ Соп@ поп 

У/ив Фе шиодисйоп о {есбпо]огу ш Фе Пееё ро\ег р!ап {есбтса| сопа1оп, опПпе 
топКопие Ваз Бесоте а {юртса| 1554е оЁ репптапепе сопйпиоиз$ топпопие оЁ Кеу рагатеет$ 
гесагез$ оЁ Фе %аН аеу1сез. ТБе аиогз сопз4ег ап аЦегпануе тефо4 Гог деегиитто Фе 
сгапКзВай гобайоп Недиепсу ап Фе зирегсВагете трт ш Фе зужет Гог таппе Фезе1 
епоше топпоппе. ТБе тефо4 1$ Базе оп ап апа|у$1$ оЁ Фе зресга| сВагасентз@с$ оЁ у1Ьго- 
асоцз$Ис $1епа| оЁ фе сот те епоше. МитепсаПу зо]уе ше ргоМет оЁ ешттайп» Фе 
еесЕ оЁ ‘еаКазе” гези те Нот Фе ПО15сгее Еоипег Тгап$Ююгт. Тре аиог$ апа[уте фе 
рипстрМе оЁ уго-асочзис тоде фа ппретеп Ше ргоседиге ог ЕЕТ ап4 Цегануе 
тейпетепЕ а]еогит Гог Недиепсу, рБазе ап атрШаае оЁ Фе теазиге4 51опа|. ТВе 1156 оЁ 
три, ошриё апа а4лазитепе рагатеег$ оР ше тодие 1$ слуеп. Те уго-асоизис тефоа Гог 
деегтштте Фе Недиепсу спагацепзис$ ор шЙааЫе таппе епеше Баз зеуега! а4уатазе$ 
оуег тадопа|, заеоезип> Фе изе оЁ 1асботаег$. Тре Нг$ 1$ фе созё оЁ шуаПаноп апа 
зпире \ау. Убго-асоизИс авегилтаноп оЁ Фе Недиепсу свагасцетзис$ оЁ фе тодие 1$ по 
а Атес( аКегпануе © ‘аспотеетз изе4 ш ргасисе Биё сап Бе изе4 шт фе зужет оЁ сопйпиои$ 
топКопие ап Фаепозте таппе Фезе] епошез. 


Статья поступила в редакцию 02.07.2012. 
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